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文章 编号 :1005-3085(2010)05-0894-07 


不 定 整数 二 次 规划 的 一 个 新 的 分 支 定 界 算法 * 


黎 健 玲 !， 马 W, E W 
(1- 广西 大 学 数学 与 信息 科学 学 院 ， 南 宁 530004; 2 海口 经 济 学 院 基础 部 ， 海 口 570203) 
摘 要 : 本 文通 过 正 交 变换 及 止 函数 的 线性 下 方 估计 得 到 不 定 整数 二 次 规划 的 可 分 离 形式 的 连续 凸 松弛 问 


题 ， 然 后 利用 Lagrangian 对 偶 技术 导出 该 凸 松弛 问题 最 优 值 的 Lagrangian 对 偶 下 界 ， 再 结合 
和 矩形 整数 对 分 ， 建 立 了 一 个 求解 不 定 整数 二 次 规划 新 的 分 支 定 界 算法 。 最 后 给 出 了 数值 试验 对 结 


果 进 行 了 比较 。 
关键 词 : 不 定 整 数 二 次 规划 ; OPH; Lagrangian 对 偶 界 ， 全 局 解 ， 分 支 定 界 
分 类 号 : AMS(2000) 90C20 中 图 分 类 号 : 0221.2 文献 标识 码 : A 
1 引言 


本 文 讨论 如 下 不 定 整 数 二 次 规划 问题 


(IQP) min f(z) =27Qz4+cla 
s.t. Ar <b, 


ZE Zn， 


其 中 Z" 是 R” 中 的 整数 向 量 集合 ，Q@ ten xn KARTE, AkEm x nF, cE R”, DE 
R”. RAURI x 是 有 界 的 


rexXxi 人 SS {rez"|0<l<zr<u), 


其 中 ! 和 尽 是 整数 向 量 。 

整数 二 次 规划 问题 在 经 济 、 金 融 、 通 信 等 部 门 有 广泛 的 应 用 03， 多 年 来 备 受 学 者 们 的 
关注 ， 并 取得 了 一 批 研究 成 果 己 9。 不 定 整 数 二 次 规划 (IQP) 是 一 个 NP 难 的 问题 ， 许 多 算 
法 是 求解 该 问题 的 特殊 情形 ， 即 0-1 二 次 规划 问题 或 二 次 背包 问题 号 21， 相 对 而 言 ， 求 解 问 
题 (IQP) 的 算法 不 是 很 多 。Erenguc 和 Benson!) 给 出 了 一 个 分 支 定 界 算法 ， 该 算法 中 下 界 的 确 
定 是 基于 正 交 变换 和 凸 包 络 技术 。Thoai 出 给 出 了 一 个 新 的 分 支 定 界 算法 ， 其 中 下 界 由 求解 凹 
规划 确定 ， 分 支 则 采用 超 和 矩形 整数 二 分 法 。 孙 小 玲 、 黎 健 玲 等 03] 讨论 了 问题 (IQP) 的 几 种 下 
界 及 相互 关系 。 
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2 A5 Lagrangian 对 偶 界 
考虑 问题 (IQP) 的 连续 松弛 
(P) min jz)=zTQz+crz 
s.t. Ax <b, 
ce XS {R"|0<l<aK<u}. 
显然 问题 (已 ) 的 最 优 值 是 问题 (IQP) 最 优 值 的 下 界 。 下 面 考虑 问题 CP.) 最 优 值 下 界 的 估计 。 
不 失 一 般 性 ， 设 和 ; < O48 = 1,---,p), à > OG = pl, n), w, wn EWN 


KEE Ht, RE ERE W = (w, Wn) WWW = En (En Anh B&H), WTQW = 
diag( 和 1,… ,和 n)。 作 正 交 变换 zz = Wy, W(P) 可 等 价 化 为 如 下 可 分 离 形 式 


(Py) min A(y) = 》 hi(yi) = 》 (čiyi + Aay?) 
i=l i=1 


s.t. Ay <b, 
l<Wy<u, 
其 中 
hi(yi) 2 CYi + NY či = Twi, i= 1,2, vn, A = AW. 

HPA <0 =1,---,p), AAIE (P,) 仍 是 非 凸 的 。 

对 于 i = 1,… ,p， 设 Qi, 记分 别 是 入 的 下 界 和 上 界 (ai, 6 的 计算 方法 参见 文献 [12])。 由 
TA <OG=1,---,p), PW Ay) (i = 1,--- ,p) SEHR. B lily) SM BR hi(y;) 在 [as Gil 
上 的 线性 下 方 估 计 ， 则 经 计算 得 

li(yi) = ey; + àa? + à:(bi + ai)lyi — ai), ye [ai, Bil, i=1, ,p. 

H Ls (ys) 取代 问题 (Py) 目标 函数 中 的 hi(yi) 得 到 如 下 可 分 离 形式 的 凸 松弛 问题 

P n 

(CP,) min h(y) = 》 lily) gz > hi(yi) 
i=1 i=p+1 
s.t. Ay < b, 
L< Wy< u. 

24 He BE n AAKRE E E (CPR Arte, FEH Lagrangian 对 偶 策 略 作 进 一 


步 分 析 。 记 w = (€7,y7,v7)7, EPE ERT, pe R, veR MEXR 知 ， 对 任意 
的 we RPH”, HE (P,) 的 对 偶 函 数 为 


Pp 
d(w) = —€Pb+uTl—vTut S min {mi(y) + (ATE + wi (v ~)) ys} 
i=l ft hs 


n 
+ in {hi(ys) + (ATE + WE w ~ u) Ty}, 
之 mia (yi) + (A E+ WE (vp)) yx} 
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WF i=1,2,---,p, hilyi) + (ATE + wT — nw)? ys E [oi p] LEM, FRA 


atin, alui) + (APE + wi W — u) ur} 


= min {hi(ai) + (APE + wp (v = u)) ai, hi (Bi) + (ATE T wi (v = u))T Bi}. 
对 于 ;=p+1 jn, hi(ys) 十 (ATE 十 wT(v — u) Ty ELZA, FRB 
min {hi(ys) + (AFE + wi (v — p)" vi} 


WATE + oP -T EP 


t 


= min {Ay + (ATE +w (v — p)” + Gui} = 
由 弱 对 偶 定 理 知 d(w) < v(P,)， 进 而 d(w) <vJQP), Ho.) RAB (.) 的 最 优 值 。 于 是 为 
得 到 最 好 的 下 界 ， 我 们 考虑 Lagrangian 对 偶 问 题 为 
(D) max d(w). 
wERT +2" 
由 文献 [12] 知 下 面 结论 成 立 。 

定理 2.1 v(CP,) =v(D)。 

上 述 定 理 表明 连续 凸 松弛 产生 的 下 界 与 Lagrangian i FRESH. (AF SHR 
大 时 解 凸 二 次 规划 (CP,) 较 费时 ， 所 以 在 下 面 给 出 的 算法 中 我 们 是 利用 解 Lagrangian 对 偶 问 
题 (DD) 确定 下 界 。 

因为 对 偶 函 数 dw) ERE PEM) eR, ATOM al i (D) 是 非 光滑 优化 问题 ， 可 用 次 梯度 
法 求解 。 


3 ”分 支 定 界 算法 
本 节 我 们 详细 给 出 求解 问题 (IQP) 的 新 算法 。 新 算法 是 一 个 分 支 定 界 算法 ， 问 题 (IQP) 最 
优 值 的 下 界 由 上 一 节 导 出 的 Lagrangian 对 偶 界 确定 ， 上 界 则 由 算法 迭代 过 程 中 产生 的 可 行 点 的 
目标 函数 值 给 出 ， 而 区 域 的 剖 分 则 使 用 超 所 形 整数 二 分 法 内 。 为 方便 起 见 ， 下 面 对 超 矩形 整数 
二 分 法 作 简 要 介绍 。 
3.1 超 矩 形 整 数 二 分 法 
WP = (a,b) C R" 是 一 个 有 整数 边界 的 超 矩 形 ， 即 
P = {xER"|a<z<b} 
= {x= (z1, an)" ER" | aj < zj < bj, j=1,2, ,n}, 
其 中 oj bj (7 =1,2, …,m) 是 整数 。 令 ;为 当前 迭代 的 分 支 指标 ， 将 书 剖 分 为 两 个 新 子 箱 
P! = {xzePl|a <r <n}, P? = {re P|n+1< r <b}, 
其 中 w= [(a; + b;)/2), t 表示 不 大 于 + 的 最 大 的 整数 。 显 然 有 


PIUP2CP PnP2=1 (P!NZ")U(P2NZ") = PNZ", 
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称 Pl，P? 为 PP 的 整数 对 分 。 

3.2 ”算法 

在 新 算法 中 ， 每 个 子 问 题 ( 子 箱 ) 对 应 于 搜索 树 的 一 个 节点 ， 通 过 解 初始 箱子 (1,w) 上 的 对 
侦 问 题 得 到 初始 下 界 。 算 法 用 一 链表 保存 子 箱 ， 进 入 第 次 迭代 时 ， 从 该 链表 选取 下 界 最 小 
的 子 箱 并 剖 分 为 两 个 新 子 箱 ， 对 每 个 新 子 箱 ， 通 过 解 对 偶 问 题 得 到 下 界 ， 若 该 下 界 小 于 问 
题 (QIP) 的 上 界 ， 则 将 该 新 子 箱 存 入 链表 中 。 利 用 剖 分 过 程 中 产生 的 可 行 点 的 目标 函数 值 更 新 
上 界 。 该 过 程 一 直 重 复 下 去 直到 链表 等 于 空 集 为 止 。 

算法 A 

步 又 0 (初始 步 ) 给 定 精度 es > 0. WHIM Q KIE A An SP = (l,u), AK 
梯度 法 解 对偶 问 题 (D)， 并 令 初 始 下 界 LB = v(D)。 车 初始 箱子 P 的 某 个 顶点 为 可 行 点 ( 记 
Aa), WAM ERUB = flat), S = {z!}， 否 则 令 初 始 上 界 UB = +0, S =. 
ST = {Pi}, k=1。 

FRID) 利用 超 和 矩形 整数 二 分 法 将 天 剖 分 为 两 个 子 箱 PL, PR, JAT HRR Pho 
令 i= 1。 i 

步骤 2( 定 界 ) 求解 对 应 于 Pi 的 对 偶 问 题 (D) 得 下 界 LB(Pi)。 若 LB(Pi) <UB-—e, W 
将 及 添 加 进 荆 中 。 利 用 剖 分 过 程 产生 的 问题 (IQP) 的 可 行 点 更 新 集合 S。 若 i =2 有 日 S Zo, 
则 令 UB = min{f(z) | zx € S}, HRH € S, 使 f(z*+1) = UB， 转 步骤 3， 否 则 令 ? = 
1 十 1， 重 复 步骤 2。 

步骤 3 (终止 性 判断 ) 将 下 中 下 界 大 于 当前 上 界 UB 的 子 箱 删除 。 如 果 工 =O, ik, ott) 
是 问题 (IQP) 的 最 优 解 ， 否 则 从 工 中 选取 具有 最 小 下 界 的 子 箱 ( 记 为 Phri)e Ok:=k4+1, 1K 
回 步骤 1。 

注 3.1 在 上 述 算 法 的 步骤 1 中 ， 我 们 沿 子 箱 最 长 边 的 中 点 经 取 整 后 对 子 箱 进行 前 分 。 

注 3.2 AE (IQP) 的 可 行 集落 在 初始 箱子 (1,w) 中 ， 所 以 经 有 限 次 超 和 矩形 对 分 后 所 得 子 
箱 的 某 些 顶 点 一 定 是 可 行 点 ， 从 而 算法 经 有 限 次 迭代 后 步骤 2 中 的 S 必 非 空 。 

因为 和 是 有 界 集合 ， 所 以 根据 整数 规划 分 支 定 界 算法 理论 ， 我 们 不 难 证 明 下 面 结论 成 立 。 

定理 3.1 算法 A 经 有 限 步 迭 代 后 可 得 到 问题 (IQP) 的 最 优 解 。 


4 ”例子 
本 节 通 过 一 个 小 例子 来 说 明 算法 A 的 实现 过 程 。 
例 
min f(x) = ~—72 — 2x22 
s.t. © + 2a. <5, 


0<zı <5, O0<22 <2, 





t= (z1, z2)” € 2Z2, 
该 问题 的 最 优 解 为 z* = (1 2)7， 最 优 值 为 f(x*) = -8。 此 问题 的 可 行 域 如 图 1 所 示 。 
第 一 次 和 迭代 
步骤 0 给 定 s = 0.001, 1 = (0,0)7, u= (5,2)7。 经 计算 得 特征 值 


Al = —1.6185, Ag = 0.6185. 
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$P = (1,w)， 解 对 偶 问 题 (D) 得 到 初始 下 界 LB = v(D) = -9.529。 在 解 对 偶 问 题 过 程 中 得 到 
HITA at = (2, 1)7， 于 是 令 5 ={(2,1)7}, 初始 上 界 UB= f(z!) =-5. ST = {Py}. 

步骤 1 箱子 的 平行 于 za 轴 的 边 较 长 (长 度 为 5)， 所 以 沿 该 边 在 7 = [248] = 2 处 应 用 整 
数 二 分 法 将 Pi 剖 分 为 下 列 两 个 子 箱 


Pi = ((0,0)",(2,2)"), P? = ((3,0)", (5,2)”), 
解 相应 的 对 偶 问 题 得 下 界 LBP) = —8.895 #1 LB(P2) = —5.178; 由 于 
LB(P}) <UB, LB(P?) <UB, 
因此 把 Pi, PRIMAT H, AEP, 从 工 中 删除 ， 于 是 工 变 为 
r = {TP} U {PI Pr} = (Pr, PP} = {((0, 0)", (2, 2)7), ((3,0)7, (5,2)7)}. 


步骤 2 在 此 过 程 中 不 需 更 新 上 界 UB， 即 UB = -5。 从 工 中 取出 下 界 最 小 的 子 箱 PH, 
£ Pz = P} = ((0,0)7, (2, 2)7). 
该 次 迭代 过 程 见 图 2 所 示 。 














dt P 的 两 条 边 一 样 长 ， 我 们 沿 平行 于 zi 轴 的 边 在 7 = 42] = 1 处 应 用 整 
数 二 分 法 把 Po 剖 分 为 下 列 两 个 子 箱 
P} = ((0,0)7,(1,2)7), P2 = ((2,0)7, (2,2)7). 
Po AAT PBR, PET BEA 
r =D\{Po} = {PP} = {((3, 0)”, (5,2)")}. 


步骤 2 PRAIA TB Le FFL B(P2) = —7.736 Al LB(P3) = —4.228; HF LB(P}) < 
UB, LB(P2)>UB, PRA A RE PI 添加 进 工 中 ， 于 是 工 变 为 


r =T|J{P2} = (Pz, Pr} = {((0,0)7, (1,2)7), (3, 0)”, (5, 2)7)}. 
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在 此 过 程 中 得 到 可 行 点 信 = (1,2)”， 所 以 集合 5 更 新 为 


令 新 上 界 


Sr? =%. KILB(P}) > UB, LB(P?)>UB, MAK PL, P? AT PAR, FRE = 0。 





$= SU {G,2)7} = {(2,1)7,0,2)7}, 


UB = min { f(x) |z € S} = f(#) = —8, 


步骤 3 AP =O, MURERE, r? = (1,2)7 即 为 问题 (IQP) 的 最 优 解 ， 最 优 值 为 fl?) 
= 一 8。 该 次 迭代 过 程 见 图 3 所 示 。 


5 ”数值 试验 





图 3: 第 二 次 迭代 





为 了 检验 算法 A 的 可 行 性 和 有 效 性 ， 我 们 使 用 Fortran 90 编写 程序 在 计算 机 Aspire-4710G- 


4A0508CIPC 上 运行 ， 


该 计算 机 的 配置 为 IntelCoreIM 双核 处 理 器 T2450 (2.0 GHz, 2MB L2 ca- 


che) 以 及 1.0GDDR2 内 存 。 测 试问 题 是 随机 生成 的 。 在 数值 试验 中 ,我们 将 本 文 的 算 
法 A( 简 记 为 Algo A) 与 文献 [各 中 的 算法 2( 简 记 为 Algo B) 作 比较 ， 数 值 结果 见 表 1， 表 中 的 
平均 CPU 时 间 和 平均 迭代 次 数 是 测试 了 10 个 问题 得 到 的 。 数 值 结果 表明 我 们 提出 的 算法 是 有 
效 和 可 行 的 ， 能 在 合理 的 时 间 内 求解 中 等 规模 的 (IQP) 问题 ， 对 同一 问题 ， 算 法 A 的 运行 时 间 


比 算法 B 少 ， 且 能 求解 的 问题 的 规模 要 比 算法 B 所 求解 的 问题 规模 大 。 


n 
5 
10 
15 
20 
25 
30 
35 


表 1: 数值 结果 
=a sme 
3 0.137 0.112 4 4 
6 0.546 0.462 8 6 
7 1.746 2.813 8 6 
13 22.25 18.17 11 15 
8 = 27.33 = 122 
8 3 33.27 5 97 
7 = 60.02 z 116 
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Abstract: With the help of the orthogonal transformation and the linear underestimation of a con- 
cave function, we first obtain a separable and continuous convex relaxation problem of the indefinite 
integer quadratic programming, and then the Lagrangian dual lower bound of the optimal value of the 
convex relaxation is derived by using the Lagrangian dual technique. Combining with integral hyper- 
rectangular bisection, a new branch and bound algorithm is presented. Finally, preliminary numerical 
results are reported. 
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